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Abstract 

It is well known that the joint probability density of the eigenvalues of Gaussian ensembles 
of random matrices may be interpreted as a Coulomb gas. We review these classical results 
for hermitian and complex random matrices, with special attention devoted to electrostatic 
analogies. We also discuss the joint probability density of the zeros of polynomials whose 
coefficients are complex Gaussian variables. This leads to a new two-dimensional solvable gas 
of interacting particles, with non-trivial interactions between particles. 

To appear in the Proceedings of the International Conference on Strongly Coupled Coulomb 
Systems, Saint-Malo, 1999 (EDP Sciences, Les Ulis). 



1. INTRODUCTION 

The study of high lying states in the spectrum of atomic nuclei led to the introduction of a 
new type of statistical mechanics H], |2|. Contrary to the traditional one, this new probabilistic 
theory does not associate a weight function to the different configurations a system may take, 
but rather introduces a statistical law directly at the level of the Hamiltonian. Depending on 
the symmetries of the system, an ensemble of Hamiltonians is defined with a proper probability 
distribution function. The statistical properties of the eigenvalues and eigenvectors of such ran- 
dom matrices (or averages of physical quantities over such distributions) may then be compared 
to a specific system, in which the average is made over some energy window. 

Motivated by its mathematical tractability, one of the original ensembles introduced is the 
so-called Gaussian ensemble of random matrices, defined as the set of matrices whose elements 
are Gaussian independent variables. The probability density of a given N x N hermitian matrix 
is defined as 

P{H)=Mexp(^-^tT{H^)y (1) 

where is a normalization constant. A very specific aspect of this ensemble is its invariance 
under rotations in Hilbert space. The form of H depends on symmetry considerations: its 
matrix elements are real for even spin systems with time reversal symmetry, complex when no 
symmetry is present and quaternion for odd spins with time reversal symmetry. For these three 
cases, the parameter f3 in Eq.(|^) takes the value 1, 2 and 4, respectively. 

During the 40 years following its introduction, random matrix theory (RMT) has proven to 
be a pervasive theory (see e.g.Q). Its applicability ranges from the description of atomic and 
molecular systems , the behavior of classically chaotic systems |Q] , the closely related motion 
of electrons in a disordered potential |5| , up to the description of the statistical properties of the 
critical zeros of the Riemann zeta function Q. In the second section of this paper our aim is to 
present a brief survey of some elementary results on the statistical properties of eigenvalues and 



eigenvectors associated to Eq.(|l|). The most characteristic feature concerning the distribution 
of eigenvalues is the existence of a strong repulsion acting among them. This repulsion may be 
seen, in an electrostatic interpretation, as a two-dimensional Coulomb interaction. For short 
distances, this Coulomb interaction is not a particular feature of the ensemble (|^) but rather 
constitutes a very basic and general fact of any quantum mechanical system. The general 
arguments leading to this result from standard perturbation theory are also presented in section 
2. 

Because of the symmetry of H, the eigenvalues lie on the real axis of the complex plane. A 
fully two-dimensional Coulomb gas arises when the basic hermitian symmetry of the Hamiltonian 
is abandoned. We illustrate this in section 3 with the example of the ensemble of complex 
matrices introduced by Ginibre Q. We discuss in some detail the nearest neighbor spacing 
distribution in the complex plane since this quantity has been less studied in two dimensions. 

In the same way the spectrum of random matrices may be associated to the equilibrium 
properties of one or two-dimensional Coulomb systems, it is possible to associate to the eigen- 
states of random matrices a two-dimensional gas of interacting particles (section 4). Although 
it contains a Coulomb repulsion, the full interaction between particles is now more complicated 
and includes many-body terms. 

This reduction of the statistical properties of eigenstates to a two-dimensional interacting 
system is done through the Bargmann or coherent state representation of quantum mechanics. 
For finite (but large) matrices like the ones considered here, this representation associates to 
each eigenstate a polynomial whose coefficients, according to RMT, are random independent 
distributed. The reduction is completed when looking at the distribution of the roots in the 
complex plane of these random polynomials. The joint probability density of the zeros may be 
interpreted as a two-dimensional gas with the properties mentioned above. 

The study of the distribution of roots of random polynomials has been an intensive field 
of research in recent years P]- pl[ . For a Gaussian distribution of coefficients the problem has 
been solved. The density ||lO| , [l2| and all higher correlation functions |13] have been computed 
analytically for arbitrary variances. The gas associated to the zeros of random polynomials 
with Gaussian coefficients constitutes therefore a further example in the restricted class of two- 
dimensional solvable models (at a certain temperature). Section 4 is a brief introduction and 
summary of some results on roots of random polynomials. 



2. RANDOM HERMITIAN MATRICES 

Consider the ensemble of NxN random hermitian matrices defined in Eq. @) . The basic problem 
to solve is to find the distribution function of the associated eigenvalues and eigenvectors. We 
denote E^, a = 1,. . . ,N the eigenvalues and {pk} a parametrization of the eigenvectors of H 
(the number of parameters pk depends on the symmetry of H). In the basis where H is diagonal, 
tvH^ = J2 ^a- Moreover, the Jacobian of the transformation is given by Q 

a 

The eigenvector's components are absent in this expression because of the rotational invariance 
of the ensemble. From these results, we obtain for the joint probability density of the eigenvalues 
(in an abuse of notation, in this paper we denote all distribution functions with the same symbol 
P) 



i<i=l V i=l / 



(3) 



j<j= 

where Z is a normalization constant. 

Interpreting the energies as the (real) positions of some fictitious particles located at Xi = Ei, 
Eq.(|^ may be interpreted as the Boltzmann factor of a classical system, P{xi, . . . ,xn) = 



Z ^ exp(— /31/(xi, . . . , Xjv)), with 

V[xx, . . .,xn) = -^log|xi - Xj\ + ^Xl^? ■ (^) 

i<j i 

The correlations among eigenvalues are therefore described by a two-dimensional Coulomb sys- 
tem made of unit charges (electrons) restricted to move over the real axis. In RMT, the temper- 
ature takes the discrete values /3 = \jkT = 1, 2, 4. The harmonic oscillator confining potential is 
interpreted as the potential produced by a uniform background of opposite charge density. This 
is therefore a one-dimensional one component plasma (IdOCP), with particles' interactions as 
in two dimensions. 

The correlations introduced by the Jacobian that are responsible for the Coulomb repulsion 
among eigenvalues are not a special feature of the particular ensemble of matrices considered 
above. In fact, for short distances compared to the mean level spacing two neighboring eigenval- 
ues of a generic quantum system repel each other as \E-^- — E^\^ . In order to see this, consider the 
energy surface of two nearby eigenstates E± as a function of a set of parameters {Ri, R2, ■ ■ ■}■ 
Two-dimensional degenerate perturbation theory states that in order to produce a degeneracy, 
(/? + 1) parameters should generically be varied. Moreover, near the degeneracy the contact 
between the two surfaces has generically the geometry of a diabolical point (a double cone) 



/3+1 



It follows from this local geometry that for short distances the probability p{s) to find two 
neighboring eigenvalues a distance s = — E^\ apart from each other is proportional to s^, 
as stated above. 

From Eqs.(|^ and (^) we see that what RMT does is to extend the simple Coulomb interaction 
which operates at short distances to all levels, located at arbitrary distances. It has moreover 



been shown |22] that the precise form of the confining potential may modify the density of 
eigenvalues, but not the local correlations between eigenvalues (at the scale of several mean 
spacings) that are therefore, in this respect, universal. These universal correlations have been 
conjectured to describe the local correlations of eigenvalues of fully chaotic systems |23, 

The partition function Z{f3) = J ■ ■ ■ J dxi . . .dxj\f exp(— /?F) of the IdOCP has been com- 
puted for any temperature Q 



N 

nr(i-/3j/2) 

27r j=i 



N 



It as also been shown that the minimum of the potential (§) (i.e., the zero temperature configu- 
ration) occurs when the particles are located at the zeros of the Hermite polynomials of degree 
N. 

At the temperatures /? = 1,2,4 required by RMT all the correlation functions of the gas, 
Rk{xi, . . . ,Xk) = {N\/{N — k)\) J . ■ ■ J dxk+i ■ ■ ■ dxj\fP{xi, . . . , xat), are known explicitly Q. For 
example, for [3 = 2 and in the large- limit they take the compact form 

i?fe(xi, . . . ,XA;) = det [sin[7r(xi - Xj)]/[7r(xi - Xj)]].^^.^-^^^^^^^ , (5) 

while the nearest-neighbor spacing distribution function is given by 

32 2 / 4 2 



P(s) = —^s exp s . (6) 

In fact, this is the 2x2 Wigner expression, that happens to be a very good approximation 
to the exact one (which takes a more complicated functional form). For applications of these 
results and a discussion of their range of applicability in chaotic and disordered systems see for 
example [0, 0]. 



3. COMPLEX RANDOM MATRICES 



Here we present an ensemble of random matrices whose (complex) eigenvalues are described by 
a fully two-dimensional Coulomb gas, as opposed to eigenvalues of random hermitian matrices 
that lie on a line. Another example of a fully two-dimensional system will be given in next 
section. 

Consider an ensemble P{H)dH of matrices whose matrix elements are independent complex 
numbers. For simplicity, we consider again a rotationally invariant ensemble 



P{H) = AAexp ( -^iT{HH 



(7) 



Again, starting from this equation the aim is to write down the joint probability density of the 
complex eigenvalues Zi,i = 1, . . . ,N . Similar steps as in the previous section leads to [0 



N 



p{zi,...,zn)=z-^ n 



z 



'j I exp 



i<j=i 



N 

E 

i=l 



(8) 



This result is formally identical to the distribution Eq.(^ for (3 = 2, with the important dif- 
ference that the eigenvalues now lie in the complex plane. Its thermodynamical interpretation 
associates to this distribution a 2dOCP (a two-dimensional gas of electrons with a uniform 
positive background). This analogy has been exploited in Coulomb systems by Jancovici |24|. 

An explicit formula for the correlations among particles with the distribution (^) can be 
worked out M 



Rk{zi,... ,Zk) = iT '"'exp 
For k = 1 and A; = 2 it reads 



El 

i=l 



det [exp(zjZj)] 



(9) 



N-l 1^12/ 



i?i(z) = iexp(-|z|2) , 
1=0 

-^2(^1, -22) = 1 - exp (-vr|zi - 22p) 



(10) 



The density of eigenvalues, that has also been studied by Girko [p^], tends to a uniform distri- 
bution on a disk of radius \z\ oc \/iV, with a Gaussian decay outside it. In the expression of R2 
the distances are measured in units such that the mean spacing is one. 

The nearest neighbor spacing distribution of the two-dimensional gas may be computed from 
the formula 



p{x)dx = 1 



1 



\Zi\>S 



d^Z2 . . . d^ZNP{zi = 0, Z2, . . . , Zn) 



(11) 



The normalization constant k is the probability to find one zero at the origin, and the notation 
is d'^z = d(Rez)d{lmz). Evaluating explicitly the r.h.s. of Eq.(ll) using Eq.(P) one finds |2£] 



/ p{x)dx = 1 - T[ \ej{s'^)exp{-s^) , 



where ej(a) = J2 01. /k\ is the truncated exponential. Deriving this expression with respect to 

k=0 

s we get 



p{s) 



-2s\[ 



r(i + i,g^) 



N-l 

E 

fc=l 



k r(fc,g^) 

r(A; + l,s2) 



1 



(12) 



where r(n, x) = /^exp(— t)t" ^dt is the incomplete gamma function. The integral of this 
function is one, as required. In order to have also the first moment set to one, one considers the 



normalized distribution Pg{s) = ap{as), with a = s p{s)ds ~ 1.14293. The smah s expansion 
of this function is pg{s) = 3Al28s^ + 0{s^). 

It is interesting to compare this distribution with the Wigner approximation in two dimen- 
sions, which gives excellent results in one dimension (cf Eq.(^). This is obtained by solving the 
approximate equation p{s)ds = [1 — p{x)dx]s'^ds. Its solution gives the normalized distribu- 
tion 

pUs) = 4 r(5/4)4 exp (-r(5/4)S4) . 
For small spacings, it behaves as p^ (s) = 2.7s^ + 0(5^) (see Fig. 1). 
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Figure 1: The nearest-neighbor spacing distribution for two-dimensional gases. Full line, Pg (Ginibre); 
dashed line, (wigner approximation); dotted line, Prp (zeros of random polynomials). 



4. RANDOM POLYNOMIALS 



As already pointed out, all ensembles of matrices previously considered are rotationally invari- 
ant. Since there is no privileged direction in Hilbert space, this induces a very simple form for 
the probability density P(ao, • • • , aat) of the coefficients {aj} of an eigenstate ip oi H expressed 
in an arbitrary basis 

P{ao,...,aN)=C 6^Ja,\^-l^ . (13) 

For convenience, we consider now matrices of dimension (A^ + 1); C is a normalization con- 
stant. We restrict moreover to the case of complex coefficients (systems without time reversal 
symmetry, j3 = 2). 

An interesting way to look at this distribution is the following. Consider the Bargmann (or 
coherent state) representation, in which the eigenstate takes the form 



j=0 



(14) 



Because we consider here finite dimensional matrices of size (A^ + 1); the usual series has been 



truncated. The polynomial ( |l4[ ) may be parametrized by its zeros 



aN 



N 



k=l 



Zk) 



It can easily be shown |27| that the statistical properties of the zeros of a polynomial are 
unchanged if the distribution (|^) is replaced by a Gaussian one 



N 



P{ao,...,aN) = Y[- — 2^^P -|«il^/(2o-, 

ZTTfJ,- 



1=0 



(15) 



The purpose now is, given the probability distribution (|15|) for the coefficients, to compute 
the distribution of the zeros of ip{z) in the complex plane. Many aspects of this problem 
for different types of polynomials (i.e., having for example some special symmetry, etc) have 
been understood in recent years |^]-|2l|. For the particular case (^) of Gaussian independent 
complex coefficients, all correlation functions have been explicitly computed [0, The 
joint probability density is easy to write down. Assuming we absorb the factors 1/a/jT in Eq.([l^) 

into the variance of the coefficients, i.e. we consider polynomials of the form ipiz) 



with the {oj} distributed according to Eq.(15). The joint probability density is then ||T 



N 
j=0 



P{zi, ...,zn) 



m 



n \Z^ 

Kj 



N 
k=0 



vr 



+ 



\zi+z2+...r 



+ 



\ziZ2+ZlZ3 + ...\ 



, \ziZ2...Zr^f 



Af+1 



(16) 



Interpreting this distribution as a Boltzmann factor, the numerator corresponds as in previous 
sections to a two-dimensional Coulomb interaction acting among the zeros. This Coulomb inter- 
action is totally generic for zeros of polynomials with complex coefficients, since it comes from 
the Jacobian of the transformation and does not depend on the particular distribution assumed 
for the coefficients. The difference with Ginibre's gas (the 2dOCP) lies in the denominator. 
The gas (|l6| ) is not confined, like the eigenvalues of random Gaussian matrices, by a uniform 
positive background inducing an harmonic well. The confinement is enforced by a much more 
complicated many body interaction between the particles. Contrary to the numerator, this term 
depends on the particular form of the distribution of the coefficients. 

In the particular case Eq.(p|) where the variances are given by cj| = the density of 

zeros in the complex plane is [Qq] 



R,{z) = I {l - g {\zf) 



l + N 



I |2 , I |2 

\z\ + \z\ g 



where g{x) = {x^ /Nl)/eNix). As for the eigenvalues of complex matrices, in the large- limit 
this density is uniform on a disk of radius oc ^/N. Outside the disk, it tends to zero much slower 
than the 2dOCP, since for large radius the tail of the density decays as l/\z\'^. 

A general expression for the two-point correlation function has been derived by Hannay [^]. 
For large N and in the uniform-density region it is given by 



-^2(^1,^2) 



f sinh^ V + v"^] cosh v — 2v sinh v 



sinh V 



with V = Tr\zi — 2;2p/2 (here again we measure the distances in units such that the mean spacing 
is one). 

A detailed comparison between the random-polynomial gas and the 2dOCP has been made 
by Forrester and Honner [ 21 1 . They found that the first moment / cP z[R2{z , 0) — 1] =— linboth 
cases, which means that a perfect screening is operating. However, the second moment differs, 
/ d^zlz^l [i?2(-2, 0) — 1] = — l/vr for the 2dOCP while it vanishes for the random polynomial gas. 



In the case of the 2dOCP, this result is known as the Stilhnger-Lovett sum rule. The physical 
interpretation of the result obtained for the zeros of random polynomials is not yet clear. 

The nearest-neighbor spacing distribution for the gas of zeros follows from Eqs.(ll) and ([T^). 
We have not been able to compute the integral. However, it may be of some practical interest to 
note that the numerical curve of p{s) for zeros of random polynomials (in the constant density 
region) is well fitted by the expression Prp{s) = 2.211s^ exp(— 0.389s^-^^^^^) (see Fig. 1). 

The validity of these results in the study of the structure of eigenstates of classically chaotic 
systems has been investigated in [12, 14, 15 1. 



5. CONCLUDING REMARKS 

This brief survey emphasizes the strong connection existing between on the one hand quantum 
statistical theories like random matrix theory or random polynomials and the physics of classical 
Coulomb systems on the other. The benefits are mutual. The thermodynamical interpretation 
serves as a guideline in the quantum mechanical studies. In the opposite direction, this connec- 
tion is a powerful tool to compute analytically the correlation functions of the classical gas, at 
least at certain temperatures. 

Many open problems lie at this interface. One example is the gas associated to the zeros 
of random polynomials, whose relevance and physical interpretation is still unclear from the 
thermodynamical point of view. This gas contains, aside a Coulomb repulsion, many body 
interactions among the particles. These interactions confine the system (they compensate the 
strong Coulomb repulsion) and also modify the correlations and the local response of the gas. 
Nothing is known about its behavior for temperatures different from (3 = 2. 

Another example concerns the fact that all the standard ensembles of random matrices or 
polynomials are related to some classical interacting gas with particles carrying the same charge 
(with eventually a compensating uniform background) . It would be interesting to study ensem- 
bles of matrices or polynomials associated to two (or multi) component plasmas, whose physical 
behavior is known to be much richer. 
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